Strained bilayer graphene: Band structure topology and Landau level spectrum 
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We show that topology of the low-energy band structure in bilayer graphene critically depends 
on mechanical deformations of the crystal which may easily develop in suspended graphene flakes. 
We describe the Lifshitz transition that takes place in strained bilayers upon splitting the parabollic 
bands at intermediate energies into several Dirac cones at the energy scale of few meV. Then, we 
show how this affects the electron Landau level spectra and the quantum Hall effect. 



PACS numbers: 73.22.Pr,62.20.-x,71.70.Di 

Electrons in bilayer graphene exhibit quite unusual 
properties: they can be viewed as 'massive chiral 
fermions' with parabolic dispersion at intermediate ener- 
gies and Berry phase 27r [1, 2], in contrast to monolayer 
graphene, where charge carriers are Berry-phase-7r quasi- 
particles with linear dispersion [3, 4]. Here, we show that 
topology of the low-energy band structure of electrons in 
bilayer graphene critically depends on mechanical defor- 
mations of the crystal. Strain determines the number of 
Dirac mini-cones in the low-energy part of the spectrum, 
below the saddle point in the electron dispersion: two 
with the Berry phases tt in a strongly strained crystal 
instead of four (three with Berry phase tt and one with 
— tt) in an unperturbed crystal [1]. These spectral fea- 
tures are tracked down to the evolution of the Landau 
levels for electrons in a magnetic field, and we predict 
their manifestation in the quantum Hall effect in strained 
bilayers. 

Bilayer graphene (BLG) consists of two honeycombs 
of carbons, ^1(2)7 ^1(2) in the bottom (top) layer and its 
band structure can be understood using the Slonczewski- 




FIG. 1: Top: Top view of the unperturbed (left) and de- 
formed (right) bilayer graphene lattice. The top and bottom 
graphene layer is shown in yellow (light) and blue (dark), re- 
spectively. Bottom: Electronic band structure of unperturbed 
bilayer graphene in the vicinity of the Brillouin zone corners 
K and K'. 



Weiss tight-binding model for graphite [5]. Within lay- 
ers, each A site is surrounded by three B sites (and 
vice versa) ^ with relative positions given by vectors ei = 

(0, l)rAB, 62 = {^,-\)rAB and 63 = {-^,-\)rAB, 
where tab is a distance between carbon sites (in equi- 
librium. Tab ^ 1-4A), and coupling 70 ~ 3eV between 
TT-orbitals of the closest carbons. The two layers are ar- 
ranged according to Bernal stacking shown in Fig. 1: 
sites A2 appear on top of Bi , with the interlayer coupling 
7i ~ 0.4eV, whereas sites Ai and B2 are positioned over 
the hexagons in the other layer and are coupled by 'skew' 
coupling 73 ~ O.SeV. All this gives rise [1] to a pair of 
'low-energy' bands touching each other in the Brillouin 
zone corners K and K' which correspond to the electron 
states located on the sites Ai and B2 , and two bands split 
by the energy ±71 which correspond to electron states 
located on the sites A2 and Bi (see Fig. 1). The degen- 
eracy of the low-energy bands is the result of cancellation 
of the contributions from the intra- and interlayer hops 
in directions ei, 62 and 63 to the transfer integral for 
the sublattice Bloch states, which is exact in K and K'. 

Perfect crystalline symmetry of graphene can be vi- 
olated when it is suspended on massive contacts: con- 
traction of the latter upon cooling can easily stretch an 
atomically thin membrane. Hence, understanding the in- 
fluence of deformations on electronic properties of BLG 
is necessary, regarding the growing interest in suspended 
graphene devices [6-10]. Below, we characterise lateral 
strain in BLG using two principal values, ^ > of the 
strain tensor u^^ = \{daUf3 + dfsUa) {a = x or y and 
u = iux^Uy) stand for displacements) and angle 9 be- 
tween its principal axes and coordinate axes in Fig. 1. 
In the case of uniaxial strain, the former correspond to 
the extension along the direction of applied tension {5) 
and contraction in the perpendicular direction {5'). Also, 
we take into account the discussed-earlier interlayer shear 
dr [11]. These deformations make the couplings 70 and 73 
dependent on the direction of the hop, which suppresses 
the above-mentioned exact cancellation in the transfer 
integral for the sublattice Bloch states and results in an 
additional term in the two-band BLG Hamiltonian de- 
scribing low-energy electrons in the vicinity of the Bril- 
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louin zone corners K = 1) and K' = —1): 



w = A3 A; Vo/3 
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^ln7o/3. 
d\n Tab ' 



As = -{S - S^)e-'''jsr]3 - ^ e^^73r/3, 

4 2 Tab 

= Pa^ + ^Pjy + Sx^A - dy^A^ 6eff, 

where 6eff plays the role of an effective magnetic field 
[12] and (j){x^ y) is determined by the condition that 
= dx'SiAQ + dy^A^^ whereas ^ = for any homoge- 
neous strain. Also, m = |^^7i/^ab7o ^ 0.035me is the 
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effective mass in the parabolic dispersion, 2^^^^ 
trons at intermediate energies 71 > e > max(m'y3, \ w\)] 
and = ^JsVab/^' The non- trivial effect of strain in 
Eq. (1) cannot be captured by the theories neglecting 
the skew coupling 73, and the effect of strain is most 
significant at low energies, |e| < max(m'U3/2, \ w\). 

To derive the effective Hamiltonian (1), we had to take 
into account that deformations modify coupling elements 
for the intralayer hops ^1(2) — ^1(2) interlayer skew 
hops Ai — B2 in directions e^, n = 1, 2, 3, 
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Here, / = (sin26>, cos 26>). These couplings enter the 
closest-neighbour tight-binding model for bilayers. 



E E (7r^4.,c.3,+H.c.) 



Here, c''" (c) are creation (annihilation) operators for elec- 
trons on the corresponding lattice sites, whereas the vec- 
tors Cji differentiate between three directions of the A—B 
hops. Note that in Eq. (2), the terms with {5 + 5') ac- 
count for 'hydrostatic' rescaling of the lattice period and 
only affect the values of m and ^3. Also, the direct A2- 
Bi interlayer coupling, 71, may be changed by shear, 
71 ^ 7i + 0{r\^)^ without any bearing on the topology 
of electron bands at low energies. 

For the electron Bloch states on the sublattice Ai, B2^ 
A2 and 5i, for the wave vectors in the vicinity of the 
Brillouin zone corners K and K' (note that strain distorts 
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the hexagonal shape of the Brillouin zone, which we also 
take into account), the four-band Hamiltonian for the 
electrons has the form 



H = 



where tt = + ipy^ v = ^vabJo/^ is the Dirac ve- 
locity in the monolayer. Following the suggestion [12] 
that in monolayers the effect of homogeneous strain ap- 
pears as constant vector potential equivalent to a small 
shift of the valley centre from the Brillouin zone corners, 
here, we employ gauge transformation of the sublattice 
spinor tj; ip ex.-p{—i^(j)) which moves the potenital part 
{dx-\-idy)(j) of Ao from the anti-diagonal part of the four- 
band Hamiltonian into the diagonal 2x2 block and add 
(subtract) -^Ao to tt {strains) in the diagonal block. Af- 
ter this, we use the Schrieffer- Wolff transformation [13] 
to project the four-band Hamiltonian onto the pair of 
low-energy bands [1] describing electron states located 
predominantly on the sublattices Ai and and, finally, 
arrive at the two-band Hamiltonian in Eq. (1). 

To judge the significance of strain, one needs to know 
the values of lattice parameters 7^0 and 7^3. Although 
the value of 773 is not known, analysis of Raman spectra 
of monolayers [14] suggest that rjo ~ —3, and we esti- 
mate that 1% of strain would lead to \w\ ~ 6meV. This 
can be further enhanced by the electron-electron interac- 
tion, which we confirm by incorporating the new strain- 
induced term in the Hamiltoniain (1) into the earlier- 
developed renormalization group theory for BLG param- 
eters [15]. The calculation, based upon the use of dynam- 
ically screened Coulomb interaction and the method of 
^/n expansion (A^ = 4 is the number of electron species 
in BLG: 2 x spin and 2 x valley), yields in the renormali- 
sation group fiow, 

dxw O.llw; dxTR-^ -0.02m"^; dxvs -0.02^3; 

where A = In ^ and e is the running energy scale. The 
electron-electron interaction enhances the strain-induced 
term stronger than other parameters, and at energies e ^ 
\w\^ where the infiuence of strain plays a dominant role 
in determining the electron spectrum, we substitute 

\w\ ^ |«;| exp(0.11 In ^ \wr^jl'\ (3) 

in the Hamiltonian (1), leading for 1% strain to an in- 
crease from 6meV to the interaction-corrected estimate, 
\w\ ^ 9meV. 

The change in topology of the low-energy dispersion for 
electrons [16] is shown in Fig. 2 for several representative 
values of homogeneous strain chosen from three charac- 
teristic regimes distinguished by shading in Fig. 3(a). 
These dispersions are plotted for the conduction band 
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FIG. 2: Calculated low-energy electronic dispersions in the conduction band of strained BLG and fan plots of Landau levels. 
Dispersion is plotted for the states near the Brillouin zone corners shown in Fig. 1 at energies |e| < lOmeV and for momentum 
in the units of p* = mvs, for several representative points in the {^w, ^w) space, as marked in Fig. 3(a). Spikes at the bottom 
of dispersion surfaces are the Dirac points characterised by Berry phases zLtt. For the Landau levels, boxed numbers mark their 
degeneracy. 



in the valley K (to be inverted in the momentum space 
to describe valley K' and flipped over for the valence 
band states, e — e). For w = 0^ electron dispersion 
undergoes Lifshitz transition [17] at the energy of the 
saddle point in the dispersion, e* = mvl/2: it splits from 
a single-connected, almost circular line into few discon- 
nected parts, each corresponding to a separate Dirac cone 
[1]. Small strain, \w\ < e*, shifts these Dirac cones across 
the momentum plane, as shown in Fig. 2 for = e* and 
w = — ze*. A stronger strain results in a collision annihi- 
lating two Dirac points, one with the Berry phase — tt and 
another with +7r, which results in a local minimum in the 
dispersion, as illustrated for w = — 5e* and w = — 3ze*. 
The other two Dirac points, each with the Berry phase 
+7r, persist to exist. In Fig. 3(a), the parametric regime 
where, in addition to a pair of well-separated Dirac cones, 
the dispersion has a local minimum, is marked by dark 
shading. Finally, much larger strain (light shading in Fig. 
3(a)) removes local minimum in the dispersion, resulting 
in even larger separation between the remaining Dirac 
cones and in a saddle point at |e| ^ \w\^ which determines 
the deformation-dependent Lifshitz transition energy in 
strained BLG. Note that all these spectral changes take 
place at a relatively low strain, ~ 1%, to compare with 



the strain of over 20% [18] required to merge Dirac points 
in monolayer graphene [19]. 

The transformation of electron dispersion by homoge- 
neous strain leads to the modification of the BLG Landau 
level (LL) spectrum. The examples of numerically cal- 
culated LLs are shown in Fig. 2 for low magnetic fields, 
B < 0.4T. Both for small and large strain, the high- 
magnetic-field end of the LL fan plot, = heB/m ^ 
max(e*, \ w\)^ is approximately described by the sequence 
e ^ ±^/n{n — l)hiJc of four- fold degenerate LLs at non- 
zero energy (n > 2) and an eight-fold degenerate LL at 
e = (n = 0, 1) [1]. In non-strained BLG at low fields, 
such that huJc{B) < rnvf, this transforms into a 16-fold 
degenerate LL at e = 0, so that the largest gap in the 
LL spectrum is between the e = and next excited LL, 
suggesting the persistence of filling factor v = ±8 in the 
quantum Hall effect (QHE) at low magnetic fields. Af- 
ter strain causes the annihilation of two out of four Dirac 
points, the e = level becomes 8-fold degenerate, and, for 
strain \w\ ^ e*, only filling factors u = +4 and u = —4 
persist in the low- field QHE in BLG: the largest energy 
gap in the LL spectra is between the 8-fold degenerate 
level at e = and next excited level, whereas the rest of 
the spectrum is quite dense. This 8-fold degeneracy is 
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FIG. 3: (a) Parametric regimes of complex w in Eq. (1), 
distinguishing between three characteristic topologies of the 
BLG spectrum in Fig. 2. (b)-(d) Activation energies for 
the QHE in strained BLG with various integer filling factors. 
For a large enough strain, (c) and (d), filling factor v — ±4 
would be the only persistent feature in the low-field QHE. 
The appearance of a local minimum in the electron disper- 
sion upon a collision of two Dirac points is manifested in (c) 
by an intermediate saturation of eact(i^) within the interval 
0.2T <B < 0.4T. 



topologically protected and it also appears in a rotation- 
ally twisted two-layer stack [20]. Figures 3(b-d) show, for 
each of the three characteristic regimes of strain in Fig. 
3(a), how the inter-LL separation (which determines the 
activation energies, eact, in the QHE) varies from high to 
low magnetic fields. At the high field end, hujc ^ 
where the LL spectrum is determined by the quadratic 
term in equation (1), eact ^ B. For the lowest 

fields, the gap between e = and the next LL scales as 
Cact oc V^, typically for the Dirac-type dispersion emerg- 
ing upon the Lifshitz transition. This behaviour is more 
pronounced for larger strain. In the regime of interme- 
diate strain corresponding to the dark area in Fig. 3(a), 
the activation energy eact of the v = A state experiences 
a very unusual intermediate saturation illustrated in Fig. 
3(c), indicating that one of the LLs gets stuck in the 
local dispersion minimum illustrated in Fig. 2. These re- 
sults are applicable locally in the case of inhomogeneous 
strain, using B^^ = B -\- ^6eff [21] which lifts the val- 
ley degeneracy of the LLs, as long as the strain varies 
smoothly on the scale longer than the effective magnetic 
length A = ^JhfeB^. 

To summarize, using homogeneous strain, one can 



spectacularly change topology of the low-energy electron 
dispersion, collide and annihilate Berry phase ±7r Dirac 
points near the corners of the Brillouin zone of bilayer 
graphene crystal. The topological changes in the disper- 
sion of electrons result in the dominance of specifically 
V — ±4 states in the QHE in BLG at low magnetic 
fields, with a characteristic behaviour of the activation 
energy as a function of a magnetic field. The latter fea- 
tures should be viewed in the context of the on-going 
experimental studies of fundamental properties of bilay- 
ers: they give a possibility to distinguish the effects of 
the deformations from spectral changes accompanying 
the earlier suggested phase transitions into ferromagnetic 
[22] or ferroelectric [23, 24] states related to the electron- 
electron interactions, which are, now, being searched for 
in suspended BLG devices [8-10]. The latter should dif- 
fer by opening a gap in the BLG spectrum. However, the 
"nematic" phase of the electron liquid [15, 25] leads to 
the same spectral changes as uniaxial strain. 
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